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Abstract 

We study a two-level atom in a double-well potential coupled to a continuum of electromagnetic 
modes (black body radiation in three dimensions at zero absolute temperature). Internal and 
external degrees of the atom couple due to recoil during emission of a photon. We provide a full 
analysis of the problem in the long wavelengths limit up to the border of the Lamb-Dicke regime, 
including a study of the internal dynamics of the atom (spontaneous emission), the tunneling 
motion, and the electric field of the emitted photon. The tunneling process itself may or may not 
decohere depending on the wavelength corresponding to the internal transition compared to the 
distance between the two wells of the external potential, as well as on the spontaneous emission 
rate compared to the tunneling frequency. Interference fringes appear in the emitted light from a 
tunneling atom, or an atom in a stationary coherent superposition of its center-of-mass motion, if 
the wavelength is comparable to the well separation, but only if the external state of the atom is 
post-selected. 
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I. INTRODUCTION 



Young's double slit experiment, in which interference is observed from light passing 
through two small slits or holes in a plate placed at a distance comparable to the wavelength 
of the light, constitutes one of the experiments at the base of quantum mechanics. Theory 
and experiment have been refined over the years to the point that the two holes have been 
replaced by two trapped atoms or ions which scatter incoming laser light 
With 



e.g. 



the advance of the coherent control of the external degrees of freedom of atoms (see 
[7, s|, 10]), the realization of atom interferometers 11, Q, Q] (see jl^] for a recent 



review), and in particular the realization of Schrodinger cat states of the center-of-mass 
coordinate of a single atom or ion [15( , it is natural to ask if interference could be observed 
in the light emitted from a single atom superposed coherently in two different positions. A 
similar question was answered to the negative in a paper by Cohen- Tannoudji et al. [3] 
for the case of scattering of light from a massive object brought into orthogonal position 
states. The physical reason is clear: interference can only arise if the probe particle can 
distinguish the two locations of the target. But then the probe particle must get entangled 
with the target. If the target is massive its two orthogonal position states remain unaltered 
during scattering, and therefore lead to vanishing overlap of the scattered probe states after 
tracing out the target. Later the scattering problem was reconsidered for li ghte r targets, 
where it was shown that interference can arise. In particular, Rohrlich et al. 17| analyzed 



the general situation of the scattering of two free particles, a probe with mass m and a 
target with mass M. Interference was predicted for the case of m ~ M, and even perfect 
visibility of interference fringes for m = M in one dimension. 

n 

Similarly, Schomerus and coworkers [18| analyzed the scattering of particles from a "quan- 
tum obstacle" , an obstacle brought into a coherent superposition of positions. An important 
difference to [l7| lies in the fact that the target was supposed to be bound in a double-well 
potential and to tunnel coherently between the two wells with tunneling frequency A. For 
the case of one dimensional scattering, they showed that the quantum obstacle leads to 
almost the same transmission resonances as two fixed obstacles, if the kinetic energy e of the 
incident particle satisfies e< A. In the opposite limit, interference can still be recovered by 
post-selecting the elastic scattering channel. 
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Spontaneous emission is not the same as scattering, and it is a priori unclear if these 
results apply to spontaneous emission alone as well. Furthermore, the properties of the 
emitted light are only a small part of the interesting physics that can arise, if internal and 
external degrees of the tunneling atom are coupled. Indeed, one might ask, if spontaneous 
emission itself (e.g. the decay rates of the excited level) changes, if the atom is brought into 
a coherent superposition of different external states. Also, what happens with the tunneling 
motion? To what extent does the emitted photon cause decoherence of the external de gree 
of freedom? Spontaneous emission from a tunneling two-level atom was considered in 19] 



for the case of the transmission of an atom through a rectangular energy barrier in one 
dimension with the atom coupled to a one dimensional mode continuum. It was shown that 
the recoil from photon emission can lift the atom over the tunnel barrier and thus increase 
significantly the transmission. 

In the present paper we examine spontaneous emission of a two-level atom trapped in a 
double-well potential, where the atom interacts with the full three dimensional continuum 
of electromagnetic modes (the interaction with a single cavity mode was treated in 20]). 
We carefully investigate the effective dynamics of all three subsystems involved: the internal 
degree of freedom of the atom, the tunneling motion, and the electric field created by the 
emitted photon in a regime where the photon wavelength is at most comparable to the well 
separation. The tunneling motion may suffer decoherence from the emission of a photon 
from the atom, but this depends, amongst other things, on the timing of the emission of the 
photon. Interference in the light emitted from the atom is very weak, but interference with 
perfect visibility of the fringes arises if the external state of the atom is post-selected in the 
energy basis. 



II. MODEL 



A. Derivation of the Hamiltonian 



We consider a trapped two- level atom (with levels \g), |e) of energy respectively) 
interacting with traveling modes of the electromagnetic field as illustrated in Fig. [TJ The 
atom is assumed to be tightly bound in the x — y plane at the equilibrium position x = 
y = and to experience a symmetric double-well potential V(z) along the z direction. The 



3 



FIG. 1: Two- level atom in a double-well potential interacting with a continuum of electromagnetic 
waves. Right panel: coordinate system used, with the atom tunneling in z-direction, and the 
atomic dipole in the x — z-plane 

Hamiltonian of this system is given by 

H = Ha + Hp + Haf-i (1) 

where Ha = H A X + if™ denotes the Hamiltonian of the trapped atom, H F is the Hamiltonian 
of the free field and Haf is the interaction Hamiltonian describing the atom-field interaction. 
Explicitly, 

a? = ^ + v(z), (2) 

flj = ^V, (3) 
H F = h^2uj k a[a k , (4) 



and, in dipole approximation, 



H AF = -J2 d - E * ( 5 ) 



k 

where d is the atomic dipole operator and 



E fe = ^e fc (a fc e tk - R + 4e- lk - R ) (6) 



the electric field operator, £ k = y ^y, eo is the permittivity of free space, V the electro- 
magnetic mode quantization volume, e k the electric field polarization vector (normalized to 
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length one), k stands for wave number k = (k x , k y , k z ) and polarization A = 1, 2 of the elec- 
tromagnetic modes with frequency = c|k| (where c is the speed of light in vacuum), and 
R = (x, y, z) for the center-of-mass position of the atom. Note that in Eqs. (E]) (JHD z is still 
an operator, with p z its conjugate momentum for the atomic center-of-mass motion along 
the z axis; M denotes the atomic mass, cr™ = \e)(e\ — \g){g\, and (at) the annihilation 
(creation) operator of mode k of the radiation field. 

In the following we will resort to the two-level approximation of the motion in the external 
potential which amounts to taking into account only the two lowest energy levels of the 
Hamiltonian H^. We denote by A the tunnel splitting, i.e. the energy spacing between the 
two lowest energy states (the symmetric |— ) and antisymmetric |+) states) of the double- 
well potential. Within this approximation, Hamiltonian (T5]) becomes 

HT = (7) 

with af- = |+)(+| — |— )( — |. We can form states that are mainly concentrated in the 
left/right wells by superposing the symmetric |— ) and antisymmetric |+) states, 



\R) 



!+>-!-> 



ex 
x 



V2 • 

The position operator z reads z = in the two-level approximation, where o\ 

+ l+)H = - \ L )( L \> and V 2 = (+\ z \-) = ( R \ Z \ R ) is the average ^-position 

of the atom localized in the right well [see Fig. [I] . 

The two-level approximation is justified if the higher vibrational energy levels are not 
populated during the spontaneous emission process. This is well satisfied, if the recoil energy 
(hw /c) 2 /2M is much smaller than the difference in energy to the next highest vibrational 
level in the external potential (Lamb-Dicke regime). For an harmonic potential, this implies 
that the wavelength A of the emitted photon is much larger than the extension of the ground 
state wave function. Here, we want to reach the regime where A ~ b to observe interference 
in the emitted light. The localization of the states \R) and \L) compared to A is measured 
by the parameter (3 = uob/c = 2irb/\. It turns out that at the border of the Lamb-Dicke 
regime in a typical double-well potential, (3 ~ 3 (i.e. A ~ 26) still leads to a reasonable 
restriction of the dynamics to the two lowest states [see Section III El for further details] . 
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By expressing the dipole operator as d = (e|d|g) cr™ = pe^a™, the interaction Hamilto- 
nian ([5]) can be cast in the form 

H AF = hJ2 9k< («^ k - R + 4e- 4k - R ) (9) 

k 

with atom-field coupling strength g^ = -pe^.e^j,/^ dipole matrix element p, and the unit 
vector Ed in the direction of the vector of dipole matrix elements, which we take without 
restriction of generality in the x — z-plane with components = (sin r/, 0, cost?) (see FigJT]). 
Furthermore, a™ = a+ + cr m , a m = \g)(e\, a+ = \e)(g\, and exp(zk.R) = cos k + ismna^ 
with k = k z b/2. We will consider in the following the situation where uj ^> A, and introduce 
the small parameter 5 = A/u <C 1. Experimentally, this is the most accessible situation 
(see Sec. Ill El) . A rotating wave approximation is then in order, which leads to the interaction 
Hamiltonian, 




Note that in the case of A ~ ujq two additional terms would have to be kept, i sin n{ak(T^a m — 
aj-crf^cr™). Equation (jTOj) makes clear that different electromagnetic waves will act in quite 
different ways: waves with sin k = will only interact with the internal degree of freedom, 
but leave the atom position untouched. Indeed, these waves do not distinguish between the 
left and the right well. Waves with sin k ^ 0, however, will couple to both internal and 
external degrees of freedom at the same time and can thus modify the tunneling behavior 
of the atom. 

It is evident from Eq. ffTUl) that the tunneling motion can in principle reduce spontaneous 
emission: The coupling constant of the second term in (111) I) changes its sign with the position 
of the atom in the double-well ("+" in the right (z > 0) well, "— " in the left well), as \R) 
and \L) are eigenstates of cx^ x with ±1 as eigenvalues. In the case of rapid tunneling motion, 
the sign of that part of the Hamiltonian is therefore reverted, so is the corresponding time 
evolution, and spontaneous emission is thus reduced. However, reverting the time evolution 
has to happen on the time scale of the dominating bath modes (i.e. 1/ujq) in order to give 
a significant effect. Based on Eq. ffTUl) . one may expect at most a reduction by a factor 2 in 
the rate of spontaneous emission for A ~ u as the first term in (1101) is independent of the 
external degree of freedom of the atom. In the limit 5 1 which we consider in this paper, 
the change of T will turn out to be very small. 
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We will describe all dynamics in the interaction picture with the free Hamiltonian 
Hq = H™ + + Hp- The corresponding time dependent field and atomic operators 
read a^it) = exp(iH t/h)akexp(—iH t/h) = a k exp(—iukt), crf£(t) = o"+ x exp(zAt), and 
= a + exp(io;oO- We thus arrive at the final form of the Hamiltonian 

H AF {t) = hJ29h(cosK (j^- u rta k o% + e^- Uk »alo^ 
k 

+i sin K(e l ^ +A - W ^a k a c *^ + e^-^^a^a^ 
B. Internal dynamics — spontaneous emission 

Let us first examine the process of spontaneous emission for the tunneling atom. We 
write a general pure state of the entire system (atom+field) as 

foK*)>= Wl n ^>, (12) 

n,<r=±,/i=3,e 

where \na/j) = |n)|(7) ex |/x)i n , |n) = rj fc \ n k) is a product state of all the field modes, n k = 
0, 1,2, . . . denotes the occupation number of mode k, and the sum over n is over all sets 
{no, ni, . . .}, |cr)ex denotes the atomic external state and |«)i n the atomic internal state . We 
start with a general initial state without any photon, but with the atom excited internally, 
and externally in an arbitrary pure state, 

1-0(0)) = c o+e (0)|0 + e) + co- e (0)|0 - e) . (13) 

Normalization imposes \c 0+e (t)\ 2 + |c — e (^) | 2 = 1- From the Schrodinger equation in the 
interaction picture, ih-^\ip(t)) = HAF(t)\ip(t)), we obtain the equations of motion for the 
relevant coefficients, 

?c -e = J2 9k e i(wo ~ Wk)t (cos k c lfc _ 9 + % sin k e~ tAt c lk+g ) , (14) 

k 

ico+e = Yl 9k e i(w °- Wfc)t (cos k c lk+g + i sin k e lAt c lk . g ) , (15) 

k 

^U-fl = 9k e- i ( w °- Wfc >* (cos k c _e i sin k e- lAt c 0+e ) , (16) 

k 

icik+a = Yl 9k e~ t(uJn)t (cos k c 0+e i sin k e iAt c - e ) , (17) 
k 
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where the dot means derivative with respect to the time t. We can formally integrate 
Eqs. ( I16|17p and insert them into Eqs. (j!4)15p . This leads to a closed system of equations for 
the coefficients Co± e . In order to avoid unnecessarily heavy notations, we focus momentarily 
on the equation for Co_ e , which can be compactly summarized as 

Co-e = G(W0, CO-e) + G(^0 ~ A, C - e ) + Gc(uo, C -e) - Gc(^0 - A, C - e ) , (18) 

£(c ,c O - e ) = ~Y,9k f dt> e^-^-^co^t') , (19) 
g c (u ,co-e) = -^^cos2k f dt'e 1 ^ xW-^co-eit'). (20) 

k ^° 

The equation for co+ e can be found by substituting A — > —A, and co- e — > co+ e in Eq. ({TBI) . 
In principle, Eq. f fT8|) contains two more terms, one given by 

I J>2 sin 2k f df e^^C^e-^^O , (21) 

and the other by an almost identical term with opposite sign and the phase factor e~* Ai ' 
replaced by e~ lAt . However, we will find that the overwhelming part of the time integrals 
comes from t' ~ t ± b/c, such that the two phases differ only by Ab/c. This quantity 
represents a tunneling speed compared to the speed of light, and has to be necessarily much 
smaller than one, as otherwise the tunneling would have to be described in relativistic 
terms. Indeed, even for very large tunneling splittings (~MHz) and well separation (~ /im, 
say), this ratio is of order 10~ 8 and thus entirely negligible. Therefore, the two additional 
terms cancel to very good approximation. 

We replace the sum over k by an integration in the limit of infinite quantization volume 
V, use polar coordinates for k, and find 

g c (u} ,c - e ) = g( + \u} ,c - e ) + g { c ~\u ,c - e ), (22) 

g( ±] {ujo,co- e ) = - aA _ 0r fc ^ / dfih(r],fi) (23) 

pt POO 

/ dt' / du^e^-^'^e^-^co^it'), (24) 
Jo Jo 



647T 2 €o^C 3 
t poo 

X 

Jo 

h(r},fi) = (sin 2 rj + 2 cos 2 r\ + /j 2 (sin 2 rj — 2 cos 2 77)) , (25) 



which is still exact, but makes Eq. (1TB]) a complicated differo-inte, 



to proceed, we resort to the Wigner Weisskopf approximation 



gral equation. In order 
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221 ] . This amounts to 



realizing that the main contribution to the integral over u will arise from a narrow interval 
around to = ujq, with a width of the order of the rate of spontaneous emission T. The 
standard Wigner-Weisskopf vacuum spontaneous emission rate reads 

r <^»- < 26 > 

which means that r(c<Jo)/^o — (4/3)a(u;od/ c ) 2 <C 1, where a ~ 1/137 is the fine-structure 
constant, and d = p/eo the dipole length (dipole matrix element divided by electron charge). 
The ratio ojq(1/c must be much smaller than one for the dipole approximation to hold. In 
our problem, the rate of spontaneous emission will be hardly modified. Thus, for all values 
of co>o, there is indeed a sharp peak of width ~ V in the integrand of the integral over u (if 
one was to perform the integration over t' first). The factor u 3 varies only slowly on that 
scale, and can therefore be pulled out of the integral. Moreover, the lower bound of the u 
integral can be extended to — oo, such that the u integral leads to a Dirac-delta function 
27c5(±fib/c — (t — £')). The slight retardation b/c corresponding to the time of travel of 
a light signal between the two wells is important here as it determines the //-interval that 
contributes, but can be neglected in co_ e after performing the integration over t', as Co- e 
will evolve on a time scale 1/T ^> b/c. We thus arrive at 

£ c (^o,c - e ) = -d(v,P)^^-co- e (t) (27) 
d( v ,/3) = - ((sin 2 v + 2 cos 2 77)^(1 - e*) + (sin 2 rj - 2 cos 2 V ) ~ 2 ' + + W - ^ 

with (3 = uob/c. The functional G(ujq, co- e ) is obtained from fl27|) by taking the limit (3 — ► 
and thus gives 

G(U , Co- e ) = -^^co_ e (t) , (28) 

as lim/3^0 d(rj, 0) = 1. Note that the dependence on u>o is both in r(a?o) and in (3. Altogether 
we have 

p 

c ±e = — ^-c 0±e , (29) 

r± = \ (tm + t(u ± A) + TMd( v , ^) - i> ± A)dfo, ^° ± A)b )) , (30) 

2 V c c / 

with the obvious solution 

C 0±e (t) = CoieWe"^ . (31) 
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It is convenient to express r± in terms of the standard Wigner-Weisskopf rate T(cu ) for a 
localized atom, and use the dimensionless parameter 5 = A/u <C 1 introduced previously. 
We then have 



For 6 = and (3 < 3 (i.e. for finite b a double-well potential with infinite barrier), we are 
immediately led back to the standard Wigner-Weisskopf rate for both initial states |0 ± e), 
r± = T(u ). Thus, as long as there is no tunneling, an arbitrary coherent superposition of 
the atom in the right and in the left well does not change the spontaneous emission at all. 
Similarly we get back F± = T(u ) for (3 = 0. For (3 > 1, d(r],j3) vanishes as 1/(3, and we 
have approximately T±/T(uo) — \ (1 + (1 ± 5) 3 ). This rate has the simple interpretation of 
arising from two independent decay channels, from | ± e) to | ± g) or | =p g). Transitions 
which do not change the external state see the same level spacing u as an atom without 
tunneling degree of freedom, whereas a flip of the external state changes the level spacing 
by ±A. Both transitions come with the corresponding Wigner-Weisskopf rate, adjusted for 
the correct overall level spacing, and the total rate is the average rate from the two decay 
channels. 

Contrary to the Wigner-Weisskopf case, the rates T± are in general complex. For the 
decay of the probabilities |co± e | 2 , only the real part of d(r],f3), 



is relevant. This function equals one for (3 = 0, depends only slightly on 77, and decays with 
some slight oscillations as a function of j3. However, for (3 > 1, 5 <^ 1 (exponentially small 
overlap of wave functions), and for S ~ 1, (3 <^ 1 (small u Q ), i.e. in the regime attainable 
within the two-level approximation made in this paper, we have that r± deviates from T(ujo) 
only very slightly and we will set r ± = T(lu ) = T in the rest of the paper. 




(32) 




(33) 
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C. External dynamics 



— decoherence of the atomic tunneling 



The quantum mechanical expectation value for the average position (z(t)) and thus the 
dynamics of the external degree of freedom follow from 

(z(t)}= b -{p-_(t)e-^ + c.c.) , (34) 

where the matrix element p^L(t) in the interaction picture is obtained from tracing out the 
internal and field degrees of freedom, 

p c + x _(t) = c 0+e (t)c* _ e (t) + K(t) (35) 

k 

The first term in Eq. (135]) can be obtained immediately from Eq. (I3TI) . In order to get K(t), 
we insert the solutions ( 1311) for co± e (t) into ( jToTfTTj) . Direct integration of the latter equations 
gives 

1 _ e (-i(w -a;fc)-r/2)t 



Cl h ±g(t) = 9k ^Co±e(0) COS K 



1 _ e (-»hTA-^)-r/2)f, 

<~^k — ^0 



1 Co ^ (0) sin K - ■•±A) + »r/2 ) • (37) 



where we have set r± = T in accordance with the previous section. It is possible to calculate 
the tunneling motion including terms of order 5, but the calculations are tedious and the 
additional information gained compared to order zero in 5 not illuminating. We therefore 
present here a simplified calculation which leads to a result valid up to corrections 0(8). 
The shifts ±A in the denominator and exponent in Eq. (|37|) have to be kept. Otherwise, 
there will be obviously no tunneling at all. Formally, the need to keep A in the denominator 
and exponent of (13T|) arises from the fact that there A has to be compared not to Uq but to 
ujq — out, which can vanish, as u>k varies from to oo. 

We restrict ourselves to real initial amplitudes Co± e (0) and proceed in a similar fashion as 
for the spontaneous emission to evaluate the sum over all modes k in Eq. (|35|) . after inserting 
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371). This leads to 



p 2 co- e (0)c o+e (0) f°° j 3 



167r 2 c 3 e ^ 



o J-l 



2/ aP,e~ rt - 2cos((c; - u )t)e- Tt l 2 + 1 
cos — 



(u - lu ) 2 + r 2 /4 



2 ./j^e 2iAt e- rt - 2e iAt cos((uj - uJ )t)e- Tt l 2 + 1 



+ sin(— ) Ai-rv a -r 

2 (a? — uj + A + i-^){uj — uj — A — % ^ 

^ 3 co- e (0)c o+e (0) p ( ~ _e- rt -2cos((a;-a;o)*)e- r */2 



167r 2 c 3 e n ./_«, I (cu-cj ) 2 + r 2 /4 

8 e 2 ^e- rt - 2e^cos((a; - ^o)t)e- rt / 2 + 1 ^ 

3 (a; — ujq + A + i^)[uj - ujq — A — i^) I 

with (3 = ujbjc and 

/ m ,-2 ,« 2^n sill ^x/'2 9 2 , /I 2/5cos/3 + (/3 2 -2)sin/3 
a(r),p) = (sin 77+2 cos 77) H 3— +(sm 77 — 2 cos 77) - + 



P J ' v ~~" '' '" V3 ' (3 3 

(39) 

We have made the same approximations as for the calculation of the rates of spontaneous 
emission, i.e. pulled out a factor Uq from the integral over u, and extended the lower bound 
of the integral to —00. In principle the cu-integral is UV divergent and would need a cut-off 
(see (23! for a discussion of experimentally relevant cut-offs). However, in order to conserve 
probability, the same approximations as for the rates of spontaneous emission need to be 
made here. Also, while there are two resonances now, uu ± A, they differ only at order 5, so 
that at lowest order in 5 it is indeed enough to pull out ujq from the integral. The remaining 
cu-integral is performed by contour integration, and we find the final result 

(z{t)) = feco-e(0)c o+e (0) e^cos(At) + ^|a(r/,/3)(l-e- r *)cos(At) 

+ (I - «v, fl) y+^Ja ( (1 + e ~ rt) sin(At) + 2 (1 " e ~ Ft) cos(At) ) } 

x(l + 0(5)), "(40) 

where we have introduced 7 = T/A. 

Let us consider a few special cases of this general result. First of all, Eq. (HOI) shows 
that for c o _e(0) = or c o + e (0) = 0, (z(t)) = for all t. This corresponds to putting the 
atom externally into one of the two energy eigenstates |±) of the uncoupled system, which 
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are symmetric with respect to z = 0. The two decay channels do not introduce a position 
bias either, and therefore the atom stays on average always at z — 0. Tunneling with full 
amplitude needs an initial preparation in the right or left well, i.e. Co_e(0)c o + e (0) = ±1/2, 
and we will therefore assume from now on Co_ e (0)co+ e (0) = 1/2. 

The limit (3 — > leads with lim^ a (Vi P) = 8/3 immediately to (z(t)) = |cos(At), 
i.e. undisturbed tunneling motion, as if the atom had no internal structure at all. The 
physical reason for this is of course that the emitted photon has in this case a wavelength 
much larger than the distance between the two wells such that it does not carry any 
information about the position of the atom, and therefore no decoherence of the tunneling 
motion arises. 



The case of f3 > 1 is more subtle. At (3 ~ 3 we have a(rj, (3) ~ 4/3 and hence 



(*(«)> - j 



+ + 778(1 



-Tt\ 



1 + 7 2 /4 



cos(At) + 



7 /4 



1 + 7 2 /4 



1 + e 



sin(At) 



(41) 



which for t ^> T , i.e. when a photon has certainly been emitted, settles down to 



(*(*)> 



7 2 A 



2 + 1 +7 2 /4 



cos(Ai) + 



7/4 



1 + 7 2 /4 



sin(At) 



with 



-Acos(At + if) 



A 



V + i 

7 2 + 4 



(42) 
(43) 

(44) 



As a consequence, after a period of initial damping, tunneling with a finite amplitude Ab/2, 
which is in general reduced compared to the full possible value 6/2, and phase shift (p persists. 

This is very much in contrast to standard decoherence scenarios of a particle tunneling 
through a potential barrier 24] . where the continued coupling to a heat-bath normally 



destroys all coherence (even though exceptions are possible in other contexts, in particular 
'or heat-baths with small cut-off frequency, which can lead to incomplete decoherence as well 



25l|). Here, the decoherence is switched off once the photon is emitted, as the center-of-mass 



coordinate of the atom does not couple directly to the electromagnetic modes. Furthermore, 
the time at which the photon is emitted, plays a crucial role. If we take 7 = V/A — > 00 
in Eq. (j4*4"I) . we find A = 1, i.e. in spite of strong dissipation and short wavelength of the 
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photon, there is no decoherence of the tunneling motion at all. The reason lies in the fact 
that the photon is emitted immediately after preparation of the atom in the right well, 
i.e. at a time, when it is not in a coherent superposition of eigenstates of its center-of-mass 
position. Thus, no coherence can get destroyed, and since after the emission of the photon 
decoherence is switched off, tunneling proceeds in the ground state with full amplitude. One 
may also see the emission of the photon as a measurement process which should, for A < b, 
project the atom either into the right or left well. But since the atom was prepared in the 
right well just before, one projects the state back into the right well, therefore keeping the 
coherence of the initial external state. 

On the other hand, if we take 7 — > in Eq. f)44p . we find that the amplitude of the 
tunneling motion in the long time limit reduces to A = 1/2. In the many runs of the 
experiment necessary to verify Eq. (j4"0|) . the time when the photon is emitted is averaged 
over many tunneling periods, such that roughly speaking in half the runs the atom is in a 
coherent superposition of eigenstates of its center-of-mass position z, half of the time it is 
in an eigenstate of z. Therefore, it is natural that on the average tunneling with half the 
full amplitude persists for t ^> 

The limits 7 — » and t —>■ 00 do not commute, which is a consequence of the factors 
exp(— Ft) in (z(t)). If we take the limit 7 — > already in Eq. (|40|) . i.e. without considering 
t — > 00 first, we get (z(t)) = | cos(At). The atom tunnels with full amplitude, i.e. A = 1, and 
shows no decoherence, as it should be, of course. Therefore, in a finite time interval starting 
with the preparation of the initial state, decoherence is most effective in an intermediate 
regime, T ~ 2 A, when a photon is likely to be emitted at the time when tunneling has 
established a coherent superposition of \R) and \L). For general 77,7 and large times, t ^> 
r -1 , i.e. after the damping has settled down, (z(t)) /(b/2) oscillates with an amplitude 



a function which we show in Fig. [2j This makes clear that also at very small (but finite 7), 
A reduces to 1/2 for sufficiently large f3. When plotting (z(t)) in a fixed time interval as in 
Fig. [21 the reduction of A is not visible yet at small values of 7. 




(45) 
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FIG. 2: (Color online) Amplitude of the tunneling motion for t S> T 1 for 7] = as function of 
and In 7. 



FIG. 3: (Color online) Tunneling motion (z(t)) in units of b/2 as function of r = At and In 7 for 
(3 = 3, 77 = with initial value (.z(O)) = 6/2. 

D. The electromagnetic field 

We now calculate the electromagnetic field in the limit of large times, t ^> Y±, such 
that the atom has emitted a photon with certainty. We will approximate again r + ~ T_ ~ 
r(u;o) = T (we are restricted to the regime uq ^> A, V). The wave function of the entire 
system has then the form 
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M°o)> 



Yl (^-ffC 00 )! 1 * - 5) + ci fc+3 (oo)|l fc + g)) 



(46) 



with 



C l fc ±5(°°) 




uj k — ojq + iT/2 



Co±e(0) COS K 



— 2 



Uk — ujq ± A + iT /2 



co T e(0) sin/« 



) 



(47) 
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(see ( 13 7p ). The first order correlation function of the electric field G^(r, r; t, t) [22J is given 
for large times by 

G«(r,r;t,t) = (V>(oo)|E<->( r , *).E<+>(r, t)^(oo)) = |I + | 2 + |r| 2 , (48) 
J± = ^c lfc± ,<0iEW(r,t)|l fc >, (49) 

with the positive frequency electric field operator E^(r, t) = ^2 k £k€kake l ^ k " r ~ ulkt \ The 
difference between I + and I - rests upon C\ k ± g . We focus for the moment on I + . The corre- 
sponding results for I~ are obtained in a completely analogous fashion. In fact, according 
to Eqs. (I47|49p one only needs to exchange c 0+e «-> c - e and replace A by —A in the final 
result to obtain I - from I + . We choose the vector r to lie in the x — z plane, convert the 
sum over k into an integral as before, and are thus led to 

I + = - I , / sin9d9 f d<p f + du^fp-k^)^^-^ 
167r 3 e c 3 J J Jo V k J 

( c o+e (0) cosK _ . c o _ e (0)sinft: \ 

X \lu - lu + iT /2 1 U-CO0 + A + iT/2J ^ ' 

= n + K , (51) 

where k.r = (uj/c)(z cos 9 + x sin 9 cos 0) . The integral splits into two parts I+, 1+ with 
denominators uj — ujo + iT/2 and uj — ujq + A + zT/2. As in the Wigner-Weisskopf theory 
of spontaneous emission, we assume that uj 3 varies little around uj = uj (respectively uj = 
ujq — A) so that we can replace uj 3 by ujq (respectively (ujo — A) 3 ) and extend the lower 
limit of integration to — oo. Only the x and z-components of 1+ and I+, denoted as 7+^, 
£ G {x, z}, give a contribution (the ^-component is identically zero), 



I+ ^opc 0+e (Q) 
c ' 5 " 167r 3 c 3 e 



-2tt 

smOdO I dcj) ft{r],6,<t>) 



o Jo 



duo — K2c '- , (52) 

f x (r),9,<f)) = [sin 7/ — sin6 l cos0(sin^cos0sin?7 + cos 6* cos 77)] , (53) 
f z (r),9,(j)) = [cos tj — cos (9(sin 9 cos (j) sin 77 -|- cos 9 cos 7])] . (54) 

Writing cos (ff COS6 1 ) as a sum of exponentials, the ^-integral can be easily evaluated 
using the contour method and leads to a Heaviside 9-function Q^t — ((2 ± 6/2) cos 9 + 
x sin 9 cos 0) /c j for the positive (negative) frequency part of the cos (|| cos #) function, 
which illustrates the fact that the electric field cannot spread faster than the speed of light. 
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In order to avoid the complications which arise from the angle dependence of the G-function, 



we will restrict ourselves to times t > -(r + 6/2) with r = \J x 2 + z 2 . We find 



b 



S7T 2 C 3 e Jo 



(55) 



xe i(o;o-ir/2)(| cos 0+f sine cos cj>-t) QQS (J-(uJ - iT/2) COS 0} 



Integration over the angle yields for the x-component 
.o;gpco+ e (0) r . 



J + = z- 
c,a: 47rc 3 e 



^ sinflcffl e <(«o-iT/2)(tcoBtf-t) cog Q_( Uq _ , r/2 ) cos ^ 



x| sin 77 J (^—(uq — iT/2) sinO^j — i sin cos 9 cos 77 Ji ^— (cjq — ^r/2) sin^ 



— sin 2 9 sin 77 



Ji(f(a;o-*r/2)siiifl) 
z(uj -iT/2) sin 9 



J 2 (-(uo-iT/2) sin 9 
V c 



(56) 



where Jq, J\, and J2 are Bessel functions. A closed form of the type of the remaining 
^-integral has been found very recently by Neves et al. [26]. Using their formula, we can 
evaluate the ^-integration analytically and find in the far-field region 



; VQpCo +e (0) i(u -iT/2)t 



47rc 3 e 
x I sin 77 
1 

COS 77 



sin R + sin R. 



ft. 



— sin 77 



sin a + sin R + sin a_ sin ft 



ft 



ft 



sin(2a + ) sinft|_ sin(2o;_) shift 



ft 



with 



ft 



ft_ 
(c^o - iT/2) 



0(1/R 



(57) 



± 



r± 



r± = \/x 2 + (z± b/2f 

X 

tan a± 



(58) 



z±b/2' 

The expression for J+ 2 is obtained from by a rotation of the coordinate system (77 — > 
77 + tt/2, a± — > a± + 7r/2), which amounts in Eq. ( IBTl) to exchanging sin 77 «-> cos 77 and 
replacing sin 2 a± — > cosa^. In the far-field, shift- ~ e lR± /2i, and the expressions can be 
further simplified, 



8llC 3 €o 



- < — — sin(77 — a + ) cosa+ + — — sin(r] — a_) cosa_ /• + 0(1/ R±) , (59) 
R+ ft 
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Similarly, we get for the second part 1+ of the integral 

r+ (go - A) 3 pc o _ e (0) i(wo _A-ir/2)t 

87rc 3 e 



x < — — sin(?7 — a+) cosa + = — sin(r/ — ot_) cosa_ f + 0(1/.R±) (60) 



with 



- (go - A - »r/2) 



(61) 



The formulas for I+ z and J+, are obtained by changing the global sign and replacing cos a± 
sina± where it appears explicitly in (1591) and (1501) . respectively . We have 



c ' x s ' a 87TC 2 e r 



(t-r/c) 



X < C O+e (0) 



e l ~ r+ sm(rj — a + ) cosa; + + e l ~ r ~ sin(?7 — a J) cosa_ 



Co- e (0) e 



Af 



sin(r/ — aq_) cos a + — e* c r ~ sin(?7 — at_) cos a 



(62) 

and again the corresponding expression for I+ z + can be found from Eq. fl62|) by just 
changing the global sign and the explicitly printed factors cosa± into sina±. We can 
summarize the results for both I + and I - as 



II 



±12 



f± _|_ J± I 2 _|_ |/± _|_ J± I 2 



w V e- r ('- ? '/ c ) 



(63) 



64vr 2 c 4 e 2 , r 2 



co±e(0) e l ~^ <5r sin(?7 — a + ) cosa + + sin(r] — a J) cos«_ 



-c OTe (0) e 



±iA(t-r-/c) 



e 



sin(r] — a+) cosa + — sin(?7 — a J) cosa_ 



+ 



c o ± e (0) e* ° 5r sin(?7 — a+) sina + + sm(r] — a_) sina;_ 



-c OTe (0) e 



±*A(l-r_/c) 



(^ =FA)ar _ 



sin(r] — a + ) sin« + — sin(?7 — a J) sina_ 

+0(A/co ,T/cu ,l/r 2 )^ 

with 5r = r + — r_. In the limit 6 — > 0, we have <5r = 0, a± = a with tana = and we 
recover the well-known result [^j (valid for t > r/c) 

G«(r, r; t, t) = -^L e" r( ^ r/c) gin 2 (r/ _ a ) (1 + 0( A /u; O) I> , 1/r 2 )) (64) 

1D7T C €q 7* 
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where rj — a is the angle between the dipole moment and the observer. 

Using 5r ~ b cos a, cosa± ~ cos a — 5a± sin a and sina± ~ sin a + 5a± cos a with 
5a± = a± — a ~ =F^ sin a in the far-field, we get for an initial delocalized state, c o ± e (0) = 1 
and c OT e(0) = 0, 

e -r(i-r/c) 



G«(r,r;M) 



16tt 2 



4^2 



ere 



sin (rj 



a) 



1 ± sm — cos a ) sm 

V2c 



Ab 

2cJ 



cos a 



(65) 



x (l + 0(A/^o,r/^ ,l/r 2 )) , 

where upper and lower sign now refer to the initial condition. 

If the atom is initially located in the right well, c o + e (0) = c o _ e (0) = l/v2 and we find 

e -T(t-r/c) 



(r,r; t,t) 



16n 2 c 4 e 2 



sin (rj 



a) 



/Ab 

1 + sin — cos a 
\2c 



cos (fl cos a) sin ( cosd 



V 2c 



— sin f A 



1 + cos \ B cos a) 



(66) 



x (l + O(A/LU ,T/LU ,l/r 2 )) . 

In both cases the interference term is proportional to sin cosa^ which scales in the non 
relativistic limit as Ab/c and is thus extremely small. 

The radiation from a classical oscillating dipole was examined recently by Bolotovskii and 
Serov 28j]. They predict interference effects in the regime where 3—1, but do not provide 
any analytical results for the visibility of the interference fringes. It appears, however, that 
a classically moving dipole has to move a distance comparable to the wavelength during the 
time 1/uq if waves from different origins are to combine in a remote location - otherwise 
the radiation pattern only follows its source adiabatically. Indeed, it is easy to show that 
at order zero in Ab/c a classical oscillating dipole does not lead to interference (where A 
means now the classical oscillation frequency of the center-of-mass motion of the dipole). 
Therefore, the absence of interference in Eq. ( 1661) at lowest order in Ab/c agrees with the 
classical result. The initially delocalized states |±) do not have a classical analogue. 

However, it turns out that in the quantum case interference with perfect visibility arises 
if the external state of the atom is post-selected in the energy basis. For example, if only 
runs of the experiments are taken into account where the atom is measured in external state 
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FIG. 4: (Color online) Polar plot of G± (r, r; t, t) (red dashed and blue dotted curves) and 
G^\r, r;t, t) = G+\r,r;t,i) + G_ (r, r;t, t) (green solid curve) for j3 = 3, rj = 0, and an ini- 
tially delocalized state co+ e (0) = 1. 

|+) before the photon is recorded, only |I + | 2 (but not |I~| 2 ) contributes to G^\r, r; t, t). We 
find in this case G+ (r, r; t, t) oc cos 2 (|cosa;) sin 2 (r/ — a) (subscript + for post-selection in 
|+)). The interference fringes for post-selection in |— ) are phase shifted, i.e., G_(r, r; t, t) oc 
sin 2 (| cosa) sin 2 (^— a), such that in the sum G^(r, r; t, t) = G+(r, r; t, t)+G_(r, r; £, t) the 
interference terms compensate and only the dipole characteristics G^'(r, r; t, t) oc sin 2 (rj — a) 
is left (see Fig. H]for rj — 0). Note that for post-selection the interference pattern becomes 
independent of A and should therefore be observable even for small A, as long as the two 
states can be reliably distinguished. 



E. Experimental perspectives 

Several requirements have to be met to observe the effects predicted in this paper. First 
of all, one needs a double-well potential with tunable well-to-well separation and barrier 
height in order to vary b and A. This has been demonstrated with optical dipole traps 
e.g. in [3, Q], and on atom chips e.g. in 0, Q]. Secondly, we considered in our model the 
same external potential for both internal states \g) and |e), and these states should be 
coupled by a dipole transition. It is by now well-known that this requirement can be met 
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for Cs, Yb, Sr, and possibly Mg and Ca atoms at certain "magical wavelengths" in optical 



traps 



30 



31 



32 



33] . Thirdly, we made the two- level approximation for the external degree 



of freedom. It turns out that in a typical double-well potential at the limit of the Lamb- 
Dicke regime, (3 ~ 3 (i.e. A ~ 2b) still leads to a reasonable restriction of the dynamics to 
the two lowest states. To show this, let us consider the usual quartic double-well potential 
V(x) = Vo(z 2 — a 2 ) 2 /a 4 where 2a is the well-to-well separation. Similar results are obtained 
for other shapes of the double-well potential. The transition probability during emission of 
a photon between energy eigenstates |n) and \m) of this potential is given by |(n|e lfcz |m)| 2 



341 ] . For Cs atoms and for the parameters Vq = 0.23 MHz and a = A/4, where A = 852.4 nm 
is the 61S1/2 — > 6P3/2 transition wavelength, this gives less than 10% transition probability 
out of the ground state doublet with a tunneling splitting A ~ 150 Hz and (3 ~ 2.93. For 
lighter atoms, such as Mg, tunnel splittings of the order of kHz can be easily achieved for the 
same (3. The theory thus works well for (3 < 3 but not anymore for larger values. Equation 
fl26|) implies that for a A in the kHz range, T ~ A for a transition in the near-infrared 
{ojq ~ 10 13 Hz), where a well-to-well separation of ~ fim leads to j3 ~ 3, which still keeps 
transition probabilities to higher vibrational states at less than about 10%. 

In any case, for optical traps the trap frequency and thus the tunnel splitting are deter- 
mined by the laser power and the focusing (or the wavelength for optical lattices), and can 
therefore be controlled independently of V such that both regimes A ^> Y and A < T should 
be achievable. The spontaneous emission rate V can be varied over a large interval by using 
a small static magnetic field to enable electrical dipole transitions between hyperfine levels 
that would otherwise be forbidden, due to admixture of small amplitudes of other hyperfine 



levels with allowed dipole transition [35], 136]. 



Cooling close to the ground state in a single well trap has been demonstrated and should 
work down to temperatures ksT < hA if A is comparable to the single well vibrational 



frequency 



301 ] . Finally, one needs to detect the tunneling motion. That should be possible 



by optical imaging, i.e. diffusion of laser light from another transition in the optical regime 
with smaller wavelength than the well separation. Another possibility might be using the 



atomic spin as a position meter 
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III. CONCLUSIONS 



A two-level atom which can tunnel between the two wells of a double-well potential allows 
for a host of interesting phenomena, which we have studied systematically in this paper in the 
regime 6 = A/u <ti 1 (tunnel frequency much smaller than the atomic transition frequency) 
and (3 = uob/c < 3. Whereas the spontaneous emission rate T of the atom is only slightly 
modified by putting the atom into a coherent symmetric superposition (ground state of the 
external double-well potential) of the two states in the right and left well, the tunneling of 
the atom and the properties of the emitted light are altered more profoundly The emission 
of a single photon can cause decoherence of the tunneling motion, but only if 1.) the photon 
wavelength is not much larger than the distance between the two potential wells, and 2.) 
spontaneous emission is not too fast, i.e. T < A. For T ^> A, the photon is emitted even 
before the atom starts its tunneling motion, i.e. before a coherent superposition of eigenstates 
of the center-of-mass coordinate of the atom is established. Hence, despite strong coupling 
to the environment, the tunneling motion does not suffer from decoherence at all in this 
regime. After the emission of the photon no more decoherence takes place, and tunneling 
will then continue with constant amplitude. For very slow spontaneous emission, r C A, 
the average amplitude of the tunneling motion reduces by a factor 2. The electric field of 
the emitted photon shows interference fringes, but their amplitude is very small unless the 
external state of the atom is post-selected in the energy basis, in which case interference 



fringes with perfect visibility arise. This is very reminiscent of the results in 18(, where 



it was predicted that interference fringes from the scattering of a particle by a quantum 
scatterer disappear in the limit where the kinetic energy e< A, but can be recovered by 
post-selecting the elastic scattering channel. The effects predicted here should in principle 
be observable with modern cold-atom technology. 
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